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ABSTRACT 

Observational studies indicate that the intergalactic medium (IGM) is highly ionized up to redshifts 
just over 6. A number of models have been developed to describe the process of reionization and the 
effects of the ionizing photons from the first luminous objects. In this paper, we study the impact of an 
X-ray background, such as high-energy photons from early quasars, on the temperature and ionization 
of the IGM prior to reionization, before the fully ionized bubbles associated with individual sources have 
overlapped. X-rays, which have large mean free paths relative to EUV photons, and their photoelectrons 
can have significant effects on the thermal and ionization balance. We find that hydrogen ionization is 
dominated by the X-ray photoionization of neutral helium and the resulting secondary electrons. Thus, 
the IGM may have been warm and weakly ionized prior to full reionization. We examine several related 
consequences, including the filtering of the baryonic Jeans mass scale, signatures in the cosmic microwave 
background, and the H~-catalyzed production of molecular hydrogen. 



1. INTRODUCTION 

Observations of the spectra of high-redshift quasars and 
galaxies shortward of Lya have revealed that the IGM is 
highly ionized up to redshifts, z ^ 6 (Fan et al. 2000; 
Hu et al. 1999; Spinrad et al. 1998). The reionization 
of the IGM subsequent to recombination at z ~ 1000 is 
thought to have been caused by increasing numbers of 
the first luminous sources, such as stars (Haiman & Loeb 
1997; Fukugita & Kawasaki 1994), quasars (Donahue & 
ShuU 1987; Valageas & Silk 1999), protogalaxies (Giroux 
& Shapiro 1996; Madau et al. 1999; Ciardi et al. 2000; 
Chiu & Ostriker 2000; Gnedin 2000), or combinations of 
these (Giroux & Shapiro 1996; Tegmark et al. 1994). 

It is unclear whether the H I and He II reionization 
epochs were coeval; it is certainly possible with a suffi- 
ciently hard ionizing source spectrum, from, e.g, quasars 
or zero-metallicity stars (Giroux & Shapiro 1996; Tumlin- 
son & Shull 2000). The present data, however, indicate 
that He II reionization occurred at z ^ 3 (Jakobsen et al. 
1994; Hogan et al. 1997; Reimers et al. 1997; Kriss et al. 
2001), and that of hydrogen before z ~ 6. The degree to 
which QSOs influence the process of reionization depends 
on their relative abundance at high redshifts. There have 
been indications that the QSO space density decreases be- 
yond a peak at z 3. This has recently been corroborated 
by the Sloan Digital Sky Survey (SDSS) observations of 
high-z QSOs (Fan et al. 2001a, b). However, QSOs may 
still be relevant to reionization, if they are powered by 
massive black holes that are postulated to form as a fixed 
universal fraction of the mass of collapsing haloes at all 
redshifts (see, e.g., Magorrian et al. (1998); Kormendy & 
Gebhardt (2001)), leading to a large population of mini- 
quasars at z > 6 (Haiman & Loeb 1998). The observed 
turnover in the QSO space density at z > 3 would then be 
true only for the brightest QSOs. 



It is also possible that some high-z QSOs have been 
missed owing to dust obscuration. The most plausible 
scenario rests on evidence that the X-ray background in- 
cludes a substantial contribution from strongly absorbed 
AGN with hard X-ray spectra (GiUi et al. 1999). Recent 
results from the Chandra satellite (Mushotzky et al. 2000; 
Giacconi et al. 2001) also suggest that many of the re- 
solved hard X-ray sources in the background could be AGN 
in dust-enshrouded galaxies or a population of quasars at 
extremely high redshift. A contrary view is provided by 
Shaver et al. (1999) who used a sample of 442 radio-loud 
quasars, unhindered by dust selection effects, to confirm 
the rapid decline in quasar space density at z > 2.5. 

In this paper, we study the effects of an X-ray back- 
ground (XRB), including the case of X-ray photons gener- 
ated by high-z quasars, on the temperature and ionization 
of the IGM prior to full reionization. X-rays have much 
larger mean free paths than extreme ultraviolet (EUV) 
photons with energies > 13.6 eV, and can permeate the 
IGM relatively uniformly. Thus, they are capable of sig- 
nificantly heating the IGM prior to the epoch when the 
bubbles that are highly ionized by EUV photons around 
individual sources overlap. Owing both to the lower pho- 
toionization cross sections of H I and He II and to the 
reduced intensity of the radiation field at X-ray, relative 
to EUV energies. X-rays alone generally do not produce 
a fully ionized IGM. However, the primary and secondary 
ionizations created by X-rays (Shull & van Steenberg 1985) 
can produce a mild increase in temperature and partial 
reionization. 

A number of authors have examined the reheating of the 
IGM in association with or preceding reionization itself 
(Gnedin & Ostriker 1997; Valageas & Silk 1999; Shapiro 
et al. 1994; Giroux & Shapiro 1996), while others have 
studied the effects of X-rays on the high-z IGM (Sub- 
rahmanyan & Cowsik 1989; CoUin-Souffrin 1991; Madau 
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& Efstathiou 1999; Haiman ct al. 2000; Oh 2001). The 
new points of interest presented in this work include: self- 
consistent solutions for the IGM temperature and ioniza- 
tion resulting purely from an XRB; the significant effects of 
secondary electrons; and the use of more accurate calcula- 
tions of photoionization cross sections and recombination 
rate coefficients. We find that the pre-reionization universe 
may be described by a model in which the first luminous 
sources and their individual EUV Stromgren spheres are 
embedded in a prc-heatcd. partially ionized IGM, rather 
than in a cold, completely neutral IGM. However, an early 
XRB does not create significant amounts of new molecu- 
lar hydrogen in the IGM, despite the increased popula- 
tion of free electrons. This is primarily due to the strong 
photo-destruction of H~ and II2 by the near IR/optical 
(henceforth IR/0) and the far-UV (UV photons in the 
Lyman- Werner bands; henceforth FUV) backgrounds, in 
the respective energy ranges 0.755 11.2 cV and 11.2 13.6 
eV, associated with any XRB generated by QSOs. We also 
find that the enhanced electron fraction in the IGM arising 
from an XRB prior to reionization may cause an overesti- 
mation of the reionization epoch as determined from the 
cosmic microwave background (CMB). 

The plan of this paper is as follows. In §2, we outline the 
IGM heating and cooling processes that we consider, and 
the models for high- 2: XRBs. In §3, we present our results 
and discuss their implications for the cosmological Jeans 
mass, the CMB, the catalysis of molecular hydrogen, 
and the inhibition of galactic outflows. We summarize our 
results in §4. 

2. X-RAY PHOTOIONIZATION MODEL 

In this section, we detail our assumptions for the XRB 
and ionizing QSO spectra, and we describe the various pro- 
cesses used to determine self-consistent IGM temperatures 
and ionization fractions of hydrogen and helium. We fo- 
cus on the heating and ionizing effects of the X-rays on the 
initially neutral background IGM, but we do not explicitly 
follow the growth of individual fully ionized regions around 
the first QSOs or star-forming galaxies. As reionization 
begins, the volume filling factor of such highly ionized re- 
gions is small. The remaining IGM is cold and neutral, 
with only a small residual electron fraction (~ 10"'') from 
the recombination epoch (Seager et al. 2000). Highly en- 
ergetic photons emitted by these sources will propagate 
much more quickly than their associated ionization fronts. 
This implies that an epoch will exist when the thermody- 
namic properties of most of the IGM may be determined 
by these X-rays. 

We model the properties of this vohime by solving for 
ionization and thermal evolution in a uniform primordial 
medium. We consider the following processes for ioniza- 
tion solutions: photoionization, coUisional ionization, case 
B radiative recombination, dielectronic recombination for 
He I, and the coupling between H and He caused by the ra- 
diation fields from the He I 24.6 eV recombination contin- 
uum and from the bound-bound transitions of He I (pho- 
ton energies at 19.8 eV, 21.2 eV, and the two-photon con- 
tinuum with an energy sum of 20.6 eV). The last of these 
processes is incorporated following the method in Oster- 
brock (1989). The exact form of the photoionization cross 
sections for H I and He II are taken from Spitzer (1978); 



we use the fit from Vcrncr et al. (1996) for He 1. 

We also account for the effects of the secondary ioniza- 
tions and excitations of H I and He I due to the energetic 
photoelectrons liberated by the X-rays (Shull & van Steen- 
berg (1985); henceforth SVS85). These introduce a further 
coupling between the ionization equilibria of H and He. 
Secondary ionization dominates over direct photoioniza- 
tion for H I in the specialized circumstance where X-rays 
are the sole source of photoionization. A typical X-ray 
photon is far more likely to be absorbed by He 1 rather 
than H I. The ejected photoelectron, however, will ionize 
many more H I atoms than He I, as H I is more abundant. 
As a result, secondary ionizations from He I photoelectrons 
and the radiation associated with He I recombination and 
excitation are the primary sources of H I ionization. As 
pointed out by SVS85, the partitioning of a primary elec- 
tron's energy, between heating the gas and secondary ion- 
izations and excitations of H 1 and He 1, is itself a fimction 
of the gas ionization fraction, x. As x decreases, and in 
particular for x < 0.1, the photoelectron deposits more of 
its energy in collisional ionizations/excitations and less in 
heat. This x-dependence is therefore important to track 
in the pre-reionization IGM, so that the role of secondary 
ionizations for H I is not underestimated. In order to be 
consistent with SVS85, who assumed that the ionization 
fractions of hydrogen and once-ionized helium were equal, 
appropriate for the interstellar environment they were con- 
sidering, we have recast the dependence in their results on 
X11+ to a direct dependence on the electron fraction, a;e 
= ne/{nu + nue), where the number density of electrons 
rie = nnii -I- nnoii + 2nHoiii- 

We compute the thermal evolution of the IGM includ- 
ing the following processes: photoelectric heating from the 
secondary electrons of H and He, as prescribed by SVS85 
(the SVS85 solution generates less heating relative to a 
prescription where 100% of the photoelectron's excess en- 
ergy goes into heating the IGM), and heating from the 
H I photoelectrons liberated by the bound-bound tran- 
sitions or the 24.6 eV recombination continuum of He I 
(here the loss of excess energy to heat is taken to be 100% 
as further ionizations of H I or He 1 are not possible). 
Cooling terms include radiative and dielectronic recombi- 
nation, thermal bremsstrahlung, Compton scattering off 
the CMB, collisional ionization and excitation, and the 
adiabatic expansion of the IGM. The values of the recom- 
bination and cooling cx)cfficients for temperatures < 10^ K 
were taken from Hummer (1994) and Hummer & Storey 
(1998), and the heating contribution from the He I two- 
photon process was calculated using the photon frequency 
distribution given in Drake et al. (1969). For the purposes 
of comparison, we will also consider cases that do not in- 
clude adiabatic cooling, which is mimicked by artificially 
holding the IGM density constant at its value at z = 10, 
so that there is no expansion cooling for z < 10. This 
would correspond to regions of the IGM that have ceased 
to participate in the Hubble expansion, but that have not 
yet begun to collapse. 

The magnitude and spectral shape of the XRB at 
high redshift is unknown, and a completely self-consistent 
model of its formation and evolution is beyond the scope 
of this paper. Although we suggest motivations for our 
extrapolations of an XRB to early epochs, the three cases 
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below merely provide a theoretical framework in which to 
study the effects of an XRB at high redshifts. We adopt 
two basic forms for the high-^; XRB, one that is more di- 
rectly tied to early QSOs, and the other, a modified version 
of the XRB given in Madau & Efstathiou (1999). In the 
first case, the specific intensity of each QSO is taken to 
have a broken power-law form, 7^ oc v~", where a = 1.8 
for hu = 13.6-300 eV (Zheng et al. 1997) and a = 0.8 for 
hu = 300 eV - 10 keV (Washburn et al. 2001). We assume 
that only those photons with energies > 1 keV permeate 
the IGM uniformly, while the less energetic photons are 
absorbed in the individual host halos. Each QSO's spec- 
trum is further attenuated by absorption by hydrogen or 
helium in the IGM. If we take the comoving number den- 
sity of QSOs to have a peak abundance of <&qso = 10^^ 
Mpc~^ at ^ = 3 and /qso^qso at earlier epochs, then the 
mean physical (not comoving) separation between QSOs 
at a given z is: 

rfQSo = (100 Mpc) /qs'/'(1 + z)-\ (1) 

We assume then that the cumulative specific intensity 
in the IGM evolves as dqso' ^i^h the evolution of the 
source density already factored in through /qso- We nor- 
malize the unattenuated mean intensity Ii, such that an 
extrapolation down to 13.6 eV, using the QSO's intrin- 
sic spectral shape, is Iq = 10"^"'^ erg cm^^ s^^ llz~^ st~^ 
at z = 3, consistent with observational constraints (see, 
e.g., Haardt & Madau (1996); Fardal et al. (1998), and 
references therein). Realistically, this underestimates the 
corresponding X-ray contribution because the radiation at 
13.6 eV is subject to some attenuation from the IGM, even 
at z = 3. 

Combining the above factors, we may parametrize the 
evolution of the IGM-filtered specific intensity for 2; > 3 
as: 

h = 3.8 X 10^24 (/ijy/300eV)-°-* /q/q e~^- 

X [(1 + z)/Af erg cm'^ s"^ Hz"! sr-\ (2) 

where we adopt /qso = lO^O-^^^-^) from the SDSS collab- 
oration (Fan et al. 2001a), and = dQso['^HO<''L/(H°) -h 
nHeOO'i/(He'') -I- nHe+'7,y(He"'")] is the IGM optical depth, 
with Uu being the appropriate photoionization cross sec- 
tion for each species. Henceforth, we refer to the above 
XRB as case S (for the Sloan survey). 

An alternate prescription for the high-z XRB is given 
in Madau & Efstathiou (1999), in which the present-day 
XRB evolves as (1 + z)"^, weighted by an exponential cut- 
off factor to account for the decreasing source density at 
high redshift. The specific intensity of this XRB is given 
by eqn. [6] of Madau & Efstathiou (1999) in units of keV 
cm^'^ s^^ keV^^ sr~^, which, over the energy range 1-10 
keV, we take to be: 

h = 7.7 (/ii//lkeV)-°-^^ exp(-/i://40keV) 

xil + zfe^p[-iz/zcf]. (3) 

We will consider two cases of this XRB: case 1 with Zc 
= 5, as in Madau & Efstathiou (1999), and case 2 without 



the exponential cutoff factor, as in CoUin-Souffrin (1991) 
where the XRB evolves cosmologically as a radiation field. 
The assumed specific intensity of each XRB increases from 
case S through cases 1 and 2. At 2: = 9, for example, their 
respective values at 2 kcV are 4.4 xlO~^^, 1.6 xlO"^**, and 
4 xlO"^'^ erg cm~^ s~^ Hz~^ sr~-^. 

The evolution of the temperature and ionization frac- 
tions is calculated as follows. We take the background cos- 
mology to be described by, fim = 0.3, il\ = 0.7, h = 0.7, 
fib/i^ = 0.019, rcMB,o = 2.728 K, and Fhc = 0.24. The 
calculation is begun at z = 12, with an initial IGM 
temperature of 20 K and residual ionization fractions of 
X11+ = 10"** and a;Hc+ = 10~^. We consider XRBs in the 
redshift range 7 < z < 12; the lower limit is roughly the 
latest epoch from observational limits (Fan c;t al. 2001b; 
Becker et al. 2001) that can be considered as being prior 
to reionization, while at z = 12 horizon effects begin to at- 
tenuate the XRB. We demonstrate the latter by equating 
the light-crossing time for half the mean physical source 
separation (given by eqn. 1) with the age of the universe 
for XRB case S, accounting for the finite time needed for 
the X-ray photons from individual sources to permeate the 
IGM. As case S is the weakest XRB considered here, with 
the lowest source density, this yields a conservative lower 
limit to the redshift at which horizon effects become im- 
portant. Furthermore, the Eddington accretion timescale 
is ~ 4.5 xlO'' yr, assuming an average QSO radiative ef- 
ficiency of 10%. Since this corresponds to z ^ 50 in our 
adopted cosmology, by 2; = 12, sources have had ~ 4 x 10^ 
yr (many times their Eddington timescale) to form and 
generate an XRB (sec, however, Haiman & Loeb (2001), 
on how this may not prove sufficient for the formation of 
the most massive black holes that likely power the bright- 
est QSOs). 

Our numerical solution for the evolving thermodynamic 
properties of the IGM solves the nonequilibrium rate equa- 
tions and energy equation for the expanding IGM [cf. 
Giroux& Shapiro (1996), eqns. (2.1)-(2.4)]. We do not fol- 
low the nonequilibrium formation and destruction of H2, 
and we do not consider detailed radiative transfer. The 
equations arc integrated with a fourth-order Runge-Kutta 
scheme with the time steps set to be less than 3% of the 
shortest relevant thermal or ionization timescale. 



3. RESULTS AND IMPLICATIONS 

In Figures 1-3, we show the temperature, and hydrogen 
and helium ionization fractions of the IGM as a function of 
redshift for XRB cases 1 and 2. Case S was found to have 

negligible effects on the thermal and ionization properties 
of the IGM, and did not alter them appreciably, over the 
redshifts that we consider, from their values in an expand- 
ing IGM without any X-rays. We therefore do not include 
case S in the figures, but emphasize that our first result 
is that there will be no significant effects on the evolving 
IGM for a conservative estimate of the pre-reionization 
XRB. 

The plots show that, depending on the choice of the 
XRB, the IGM temperature ranges from ^ 100 K to ~ 
10^ K, with a:H+ varying from 0.1% to ^ 20%. For each 
case, the solutions including and excluding adiabatic cool- 
ing are almost identical, and are practically overlaid in 
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the figures. As may be expected, the stronger the XRB, 
the greater are the heating and ionization of the IGM. 
Though we do not show it here, an XRB whose magni- 
tude hes between those of cases 1 and 2 generates IGM 
temperatures and ionization fractions that are bracketed 
by the two cases shown in the figures. 
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Fig. 1. — The IGM temperature is shown as a function of redshlft 
for the two XRBs described in the text, each with (dashed lines) and 
without (solid lines) adiabatic cooling from cosmological expansion. 
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Fig. 2. — The IGM ionization fraction of hydrogen is shown as 
a function of rcdshift for the two XRBs described in the text, each 
with (dashed lines) and without (solid lines) adiabatic cooling from 
cosmological expansion. 

Note that for case 2, Tigm is slightly higher when adi- 
abatic cooling is included rather than excluded. This 
counter-intuitive result arises from our choice to mimic 
the cases without adiabatic cooling by stopping the cos- 
mological evolution of the IGM density at 2; = 10. For 
case 2, by z = 10, cooling is dominated by line excita- 
tion cooling, which is proportional to the square of the 
IGM density, and not adiabatic cooling. Consequently, for 
the cases including adiabatic cooling, where the IGM den- 
sity continues to decrease at z < 10, the cooling decreases 
more rapidly than it does in the absence of adiabatic cool- 
ing, leading to marginally larger values of Tigm- Note also 
that when x^, which roughly tracks Xii+, exceeds about 



10% (case 2), secondary ionization becomes less impor- 
tant for H I, and a;He+ begins to exceed Xii+ . We find that 
the He III fraction Xiie+2 was typically smaller by a factor 
of 100 than Xhs+j and we do not display it in the figures. 
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Fig. 3. — The IGM Ionization fraction of Hc+ Is shown as a 
function of rcdshift for the two XRBs described in the text, each 
with (dashed lines) and without (solid lines) adiabatic cooling from 
cosmological expansion. 

We note here some consequences of our adopted pho- 
toionization cross sections and of consistently including 
the secondary effects of photoelectrons. The He I cross- 
section is taken from the fit to experimental data provided 
by Verner et al. (1996). This differs dramatically in the 
keV energy range from the two-powcr-law fit from Oster- 
brock (1989) used by Shapiro & Kang (1987) and Abel 
et al. (1997), which is only valid for photon energies of a 
factor of a few above the threshold value of 24.6 eV. This 
power-law extrapolation exceeds the Verner et al. (1996) 
He I cross-section by a factor of 40 at 1 keV, and by a 
factor of 400 at 10 keV. It is therefore important to use 
the latter fit for He I when considering the effects of a 
pure XRB. Furthermore, the ratio of the H I/He I pho- 
toionization cross sections is very small in the X-ray range, 
about 3% for photon energies of 1-10 keV. Thus, in the 
pre-reionization universe, when the IGM is predominantly 
neutral and can be permeated only by X-rays, the ion- 
ization balance is driven primarily by the photoionization 
of He I and its associated photoelectrons (and not those 
from H I) , each of which then ionize about a dozen hydro- 
gen atoms (SVS85). Since the primordial ratio of helium to 
hydrogen is about 8%, we may anticipate that a;He+ ~ ^11+ 
in such a scenario, and that is what we find. 

For the purposes of the discussion below, we now take 
case 2 with adiabatic cooling to be our fiducial case, as it 
provides an upper bound to the magnitude of the results 
in Figures 1-3, or equivalently, it represents the "maxi- 
mal effect" of a high-z XRB. In this case, H I ionization 
is dominated by secondary ionizations from the He I pho- 
toelectrons, as expected, followed by radiation from the 
bound-bound transitions of He I, while He I is most af- 
fected by direct photoionization and secondary ionizations 
from its own photoelectrons. The major sources of IGM 
heating arc, in decreasing order, direct photoelectric heat- 
ing of He I, H I and He II, and heating from the bound- 
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bound transitions of He I. At the earliest epochs that we 
consider, adiabatic expansion and Compton cooling were 
the dominant cooling processes. For Case 2, by z ~ 11, 
the IGM had experienced sufficient heating and ionization 
that line excitation cooling was the most significant source 
of cooling. For Case 1, adiabatic expansion remained the 
dominant cooling term. 

3.1. Jeans Mass Filtering 

We now discuss the implications of such a heated, par- 
tially ionized IGM, beginning with the evolution of over- 
densities in a warm IGM. At any redshift, dark matter 
perturbations with virial temperatures Tvir(2;) will decou- 
ple from universal expansion on the associated mass scales 
and begin to collapse. The baryons will also fall into these 
dark matter potential wells if they have little pressure as- 
sociated with them. If the gas in the IGM is heated by 
luminous sources, and if Tjgm exceeds Tvir, the baryons 
will resist collapse on the mass scale set by Tigm, lead- 
ing to a filtering of the mass scale that can collapse and 
virialize. This feedback from a heated IGM has been ex- 
plored for a number of reheating scenarios (Tegmark et al. 
1997; Shapiro et al. 1994; Valageas & Silk 1999). As we 
do not account for gas clumping in our analysis, we follow 
Shapiro et al. (1994) in evaluating the magnitude of the 
filtering for the cosmological Jeans mass. Assuming that 
the temperature of the baryons and the CMB are coupled 
until z ~ 150, and that the IGM cools adiabatically after 
that, the Jeans mass in baryons is: 



Mj w 2.2 X 10^ 



h {lJ,ftm) 



1.5 



1 + Z 

10 



1.5 



M©. (4) 



If, however, the IGM is heated to a temperature Tigm, 
the Jeans mass will scale as: 



M, 



1.3 X 10^ 



Oh 



h (/iOn 



U-5 



IGM 



CMB 



1.5 



Mq, (5) 



where Tcmb = Icmb,o (1 + z)- For the purposes of com- 
parison, we work with the fiducial case defined above. At 
z = 9, Tigm ^ 1-3 xlO"* K, so that Mj - 3 x 10** Mq, 
as opposed to Mj ~ 520 Mq from eqn. 4. Thus, the 
Jeans mass in a pre-heated IGM exceeds by several or- 
ders of magnitude its corresponding value in the case of 
an IGM that cools adiabatically after decoupling from the 
CMB. Furthermore, the filtered Jeans mass increases with 
decreasing z, unlike the canonical Jeans mass. 

3.2. The CMB 

A partially ionized IGM can afi'ect the CMB through 
Thomson scattering of the CMB photons off the free 
electrons in the IGM. There are at least two poten- 
tially observable effects. First, the optical depth to elec- 
tron scattering, Te, as measured from the primary CMB 
anisotropics, will have some contribution from the par- 
tially ionized IGM preceding complete reionization, and 
may no longer be a strong constraint on the epoch of reion- 
ization. The expression for Te as a function of redshift in a 
A-cosmology is given by (see, e.g., Tegmark & Silk (1995)): 



Te{z) ~ 0.057 Ob/l X 

dz' 



J 

Jo 



{l + zTXe{z') 



V^A + (1 + Z'^il - Oa + i^mZ' 



-- (6) 



This assumes that all the baryons are in a homogeneous 
IGM at all times, which may not be a fair approximation 
at late epochs. As a simple calculation, let us examine two 
possibilities: one, that reionization occurred at z = 6, and 
that Xii+ is given by our standard XRB for 7 < 2: < 12, 
with a post-recombination residual ionization fraction of 
10^^ for 12 < z < 1000. The values of x-j^+ are linearly ex- 
trapolated between redshifts of 6 and 7. The second case 
assumes that there is no high- 2; XRB, and that xii+ = 
10~* for 7 < z < 1000. Figure 4 displays Tc as a function 
of redshift for these two scenarios. In the former case, the 
integrated optical depth to recombination is Te ~ 0.047, 
of which 0.032, or about 68%, is from the reionized IGM, 
and 0.009 comes from the partially ionized IGM, roughly 
a 20% effect. Though this seems a small correction, it is 
a level of accuracy that can in principle be probed by fu- 
ture CMB experiments such as Planck (Eisenstein et al. 
1999). The additional contribution from the partially ion- 
ized IGM could cause an overestimation of the reionization 
epoch, Zr, if the cumulative optical depth were attributed 
solely to the IGM after complete reionization. Unfortu- 
nately, the correction to Te from this XRB has roughly 
the same value as that from the IGM's residual electron 
fraction at 12 < z < 1000 (re ~ 0.006), due to the late 
reionization model and a low total Tg. 
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Fig. 4. — The optical depth to electron scattering as a function of 
redshift. Sohd hne assumes reionization of the IGM at z = 6, partial 
ionization for 7 < z < 12 for XRB case 2 with adiabatic cooling, 
and the primordial electron freeze-in abundance at higher redshifts. 
Dashed line represents the same scenario, but with no XRB. 



Since the theoretical uncertainty in the post- 
recombination a;c is of order 10% at present, all that 
is required to distinguish between the "with XRB" case 
in Figure 4 and a "No XRB" case with somewhat ear- 
lier full reionization is an independent measurement of 
Zj.. While Planck could in principle map out the low- 
multipole (l) polarization peak caused by reionization in 
the CMB angular power spectrum to significant accu- 
racy, for late reionization this feature is of suSiciently 
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low magnitude and at large enough scales {I < 10) that 
foregrounds are likely to be a complication. Ultimately, a 
spectroscopic determination of the reionization epoch via 
the Gunn-Peterson effect from the SDSS or from the Next 
Generation Space Telescope, in conjunction with future 
CMB data, may provide the best probe of an IGM par- 
tially ionized by X-rays prior to full reionization. The first 
of these observations is a real possibility in the near future 
from efforts that are well underway with data on SDSS 
QSOs at redshifts > 5.8 (Fan ct al. 2001b: Becker ct al. 
2001). Otherwise, for all practical purposes, it may be a 
challenge to separate these two different contributions to 
Te from the epochs preceding reionization. 

Although a partially ionized IGM may smear a direct 
relation between and as determined from the CMB, 
the majority of CMB photons that are Thomson-scattered 
at late times are still scattered by a fully ionized IGM. To 
sec this, consider that the probability that a CMB photon 
scattered off an IGM electron after rcdshift z along the line 
of sight is [1 - exp(-Te(z))]. The derivative of this, called 
the visibility function (see, e.g., Tegmark et al. (1994)), is 
the relative probability of the last scattering redshift of a 
CMB photon. We show the visibility fimction in Figure 5, 
for the two scenarios in Figure 4. We see that, despite the 
contribution to Te from an IGM that is partially ionized by 
X-rays, the surface of last scattering after recombination 
for the CMB photons that we see today is still strongly bi- 
ased towards the epochs following c;onipk;te reionization. 
Though we do not display it here, we note that the peak 
in the visibility function corresponding to recombination 
is much larger than that associated with reionization (see, 
e.g., Zaldarriaga (1997)), especially when the latter occurs 
at low redshifts. 




0.0 5.0 10.0 15.0 



z 

Fig. 5. The probability distribution of a CMB photon being 
Thomson-scattered by an IGM electron as a function of redshift, 
with the same assumptions as in Figure 4. 

A second signature involves the correlations between 
ionized regions in a patchy reionization model, which ap- 
pear in the multipole moments of the CMB anisotropics, 
typically at very small angular scales {I > 5000) (Knox 
et al. 1998). Since these correlations are usually predicted 
assuming a neutral background IGM, one may pose the 
question of whether such correlations may be significantly 
erased in the presence of a partially ionized IGM. As it 



happens, this is not a concern since the dominant con- 
tribution to the above correlations come from patches of 
full ionization, which are themselves highly correlated with 
high-density regions. Thus, the relevant factor is the cor- 
relation of the underlying matter density field, rather than 
the sizes of ionized (x = 1) regions or the ionization con- 
trast between x = 1 regions and the background IGM. 
The relative contributions from these different effects is 
best seen in Valageas et al. (2001). 

3.3. Molecular Hydrogen 

We now turn to the implications of an X-ray heated, 
partially ionized IGM for the production of molecular hy- 
drogen, II2, which is an important coolant for primordial 
gas in virializing structures. A prerequisite for this is that 
the baryons are able to collapse into virializing dark mat- 
ter halos, which may be hindered by a heated IGM, as 
discussed above. We emphasize that we focus here on H2 
creation and destruction in the IGM, rather than within 
dense collapsed structures. We include the following pro- 
cesses in an equilibrium calculation of the abundance of H2 
relative to hydrogen, in which we consider only the H" cre- 
ation channel and do not include formation via (Lepp 
& Shull 1984; Abel et al. 1997). Unless stated otherwise, 
rates for individual reactions are taken from Donahue & 
Shull (1991). H~ is formed through radiative attachment 
of H I (Galli & Palla 1998), and is destroyed by associative 
detachment by H L mutual neutralization with protons, 
coUisional detachment by electrons (Abel et al. 1997), and 
photodetachment by photons of energy 0.755 to 13.6 eV. 
We assume that H2 forms by associative detachment of H I 
and H", and is destroyed by collisions with electrons and 
with H I, by charge exchange with protons, and by pho- 
todissociation in the Lyman and Werner bands by photons 
of energies 11.2-13.6 eV. 

As we do not specifically account for radiation below ~ 
1 keV, we extrapolate the flux of each XRB at 2 keV down 
to 2500 A (4.97 eV) using an effective X-ray spectral index 
of {ceo,x) = - 1.38, with a power-law for photon ener- 
gies below about 5 eV. We then compute the photodissoci- 
ation rate of directly from the tabulated cross-sections 
of Wishart (1979) using the fit given in Tegmark et al. 
(1997), and of H2 as described in Donahue & Shull (1991), 
for the respective photon energy ranges stated above. For 
both and H2, particularly the latter, the photodissocia- 
tion rates exceed the other destruction channels by several 
orders of magnitude. As we have performed an equilibrium 
calculation and have not included the post-recombination 
freeze-in abundance of H2 relative to hydrogen of ~ 10~^ 
in the IGM, our results should be taken as representing 
the effect of an XRB for newly created or destroyed II2 . 

We find that, for XRB cases 1 and 2, n(Il2)/nH never 
exceeds about 3 x 10"^'*, and has little evolution over the 
redshifts that we consider; this falls well short of the crit- 
ical value of ^ 10^^ to stimulate gas cooling in high-z 
structures. Although the presence of X-rays boosts the 
free electron fraction, at any given redshift the stronger 
XRB creates less new H2 , owing to the combination of var- 
ious dissociation processes that are enhanced in a heated 
IGM, and to the strong role played by photodissociation 
in the destruction of both H~ and H2. We find that the 
timescales associated with the formation and photode- 
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Table 1 

Molecular Hydrogen Production 



Background n(H )/nH n(H2)/nH 



XRB 


1.1 


xlO 


-7 


1.6 


xlO" 


5 


XRB + IR/0 


3.3 


xlO 


-12 


4.8 


xlO" 


10 


XRB + FUV 


1.1 


xlO 


-7 


4.2 


xlO" 


10 


XRB + FUV + m/0 


3.3 


XlO 


-12 


1.2 


xKJ- 


14 



Note. — The relative roles played by X-ray, FUV (Lyman- Werner bands) and IR/optical backgrounds (sec text for the individual definitions 
of photon energy ranges) in an equilibrium calculation of n(H~) and n(H2) for XRB case 2 with adiabatic cooling at ^ = 9. 



struction of and H2 arc always much less than the 
Hubble time at each epoch, so that our equilibrium calcu- 
lations here indicate real changes to the primordial freeze- 
in abundance of H2. 

To illustrate the relative roles played by an XRB and the 
associated radiation fields below 13.6 eV, we display in Ta- 
ble 1 the equilibrium abundances of H~ and H2 calculated 
by successively including radiation from two backgrounds, 
defined as follows. The FUV includes the photons rele- 
vant for H2 photodissociation in the Lyman- Werner bands 
(11.2-13.6 eV). The H~ is dissociated by photons hav- 
ing energies of 0.755-13.6 eV, in the framework of this 
paper. For the purposes of Table 1, however, the IR/0 
(infrared/optical) term refers to the inclusion of all sub- 
13.6 eV photons except those in the LW bands. Although 
this energy range (0.755-11.2 eV) corresponds to some UV 
photons as well, most of the photodissociation occurs at 
near IR and optical wavelengths, given the nature of the 
cross section and our adopted power-law form for the 
background intensity at these energies. 

We see that the minimum IR/ O and FUV backgrounds 
associated with any putative XRB from high-z quasars 
more than compensate for any positive feedback from the 
XRB in the form of an increased electron fraction in the 
IGM (see, however, Ricotti et al. (2001) on the positive 
feedback for H2 formation in the vicinity of individual ion- 
ization fronts generated by hard stellar spectra) . The neg- 
ative feedback from FUV radiation has already been noted 
in several works; we point out here the additional feed- 
back from IR/0 photons. The importance of the IR/0 
photons in this paper relative to the results of previous 
works, e.g., Haiman et al. (2000), can be traced to at 
least two factors. Firstly, our equilibrium calculation is 
performed for the relatively low-density conditions in the 
IGM at 2; = 9, rather than within highly overdense col- 
lapsed halos, so that some variance can be attributed to 
the differing densities of the studied environments. Sec- 
ondly, Table 1 represents the effects of the strongest XRB 
considered here, which has a correspondingly high IR/0 
associated background. An XRB that was extrapolated 
from the observed EUV specific intensity at z = 3, such 
as case S in this work, would have a related IR/0 back- 
ground that has a negligible impact on the net H~, and 
hence H2, abundances. Note also that a combination of 



these two factors, IGM density and radiation intensity at 
the H" photodissociation threshold, implies that only a 
relatively modest IR/0 background is required at z < 12 
to destroy H" . A radiation field of far greater energy den- 
sity is needed to accomplish this in the early universe; this 
leads to the well-known result that the CMB photons are 
believed to have photodissociated H" at ^ > 100. 

Our estimate of the FUV background is not necessar- 
ily the minimal value, due to the opacity of the IGM in 
the Lyman- Werner bands which we have not accounted 
for here. The table shows, however, that even an order 
of magnitude reduction in the FUV background alone will 
not increase the fractional abundance of H2 above lO^*'. 
As we have neglected contributions from stellar radiation, 
we have in fact underestimated the effects of IR/0 and 
FUV photons at the epochs we consider. In summary, we 
find that an XRB does not produce significant amounts of 
new H2 in the IGM, unless the associated FUV and IR/0 
backgrounds were somehow strongly attenuated. 

3.4. Inhibition of Outflows 

Finally, it may be possible for a heated IGM with suf- 
ficient pressure to inhibit outflows from star-forming pro- 
togalaxies. Mass loss from early objects through galactic 
winds or evaporation is often invoked to explain the ubiq- 
uitous presence of metals in the Lya forest clouds at z ~ 3; 
a pre-heated IGM could, however, hinder such outflows of 
metal- enriched gas from host galaxies at sufficiently high 
redshifts (prior to reionization) . An estimate of the ef- 
fect of a warm IGM in inhibiting galactic winds may be 
made as follows. Setting the ram pressure of the wind 
PwWw> equal to Figm = (^h + riue + n^) kTjGM, and us- 
ing Mw = 47rrgj PwVw, where Mw and rgt are respectively 
the mass outflow rate of the wind and the wind stalling 
radius in the IGM where pressure equilibrium is reached, 
we have: 

r,t ^ (0.5 Mpc) X 




for the IGM density at z = 9. Note that if the IGM had 
cooled adiabatically from z ~ 150, then Tjgm would be a 
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few Kelvin, and rgt 33 Mpc. This neglects the added 
contribution that any ionizing photons escaping from the 
protogalaxy would make to the local ambient pressure. 

Eqn. 7 approximately indicates the maximum degree, 
case 2 being the strongest XRB considered here, to which 
an early XRB may hinder outflows from starforming galax- 
ies at ^; = 9. Note, however, that at lower redshifts, e.g., 
2; ~ 3, Vst can be significantly larger, particularly for the 
mass outflow rates of a few to ten solar masses per year 
that arc typical of starbursting galaxies. Thus, the obser- 
vations of a trace metallicity in the IGM at 2; ~ 3 does 
not necessarily rule out an early XRB from QSOs prior 
to rcionization, if the epoch of metal ejection associated 
with stellar activity in high-z galaxies succeeds that of 
pre-heating from an XRB by a sufficiently long period. 
This becomes a particularly important constraint, how- 
ever, for a model that posits an early XRB from stars, 
or one that explores the heating effects of X-rays for the 
post-reionization IGM (Madau & Efstathiou 1999). 

4. CONCLUSIONS 

We have examined the effects of various high-z XRBs, 

including a case which mimics X-ray photons from early 
bright QSOs, on the temperature and ionization of the 
IGM before rcionization is complete. We have found that 
individual luminous sources with their associated H II/He 
III ionized regions may be embedded in a warm, partially 
ionized IGM, rather than in a cold neutral IGM, during 
the epochs before individual Stromgren spheres have over- 
lapped. 

The heating and ionizing effects from the XRBs were 
determined self-consistently, including the associated gen- 
eration of photoelectrons which act as secondary sources, 
resulting in an enhanced population of free electrons. We 
have investigated two XRB cases that are related to the 
measured present-day soft XRB, and find that they can 
heat the IGM to temperatures between a hundred to 10^ 
K prior to rcionization, and result in ionization levels be- 
tween fractions of a percent to about 20%. The third 
XRB, assumed to be generated by early QSOs whose space 
density evolves as indicated by the recent observations of 
bright high-z quasars from the SDSS collaboration, did 
not appreciably alter the IGM temperature and ionization 
from their post-recombination values at the redshifts we 
considered. 

We have examined the implications of our results for 
several phenomena of interest to cosmology, including the 
Jeans mass, the CMB, and the production of H2. In par- 
ticular, the pre-heated IGM raises the Jeans mass signif- 
icantly at the redshifts of interest to the growth of struc- 
ture, leading to a filtering of the baryonic mass scales 
that can collapse into virializing dark matter halos. We 
also flnd that the XRB-enhanced electron fraction in the 
IGM prior to rcionization increases the total optical depth 
to electron scattering. This could lead to an overestima- 
tion of the rcionization epoch if the cumulative To as mea- 
sured by CMB experiments were attributed solely to the 
fully ionized IGM. However, the addition to from the 
post-recombination IGM, with a frcczc-in ^ 10~^ for 
z < 1000, is roughly comparable to that from an IGM 



which is partially ionized by our standard XRB. It may be 
a challenge to observationally distinguish these disparate 
contributions to Te in late rcionization scenarios, unless 
there is an independent measurement of the rcionization 
epoch, e.g. spcctroscopically through the Gunn-Pcterson 
effect from the SDSS or from future space telescopes. Such 
a constraint seems particularly promising from ongoing 
analyses that are well underway with SDSS data on QSOs 
at z > 5.8 (Fan et al. 2001b; Becker et al. 2001). In con- 
junc;tion with forthcoming CMB data, it could provide the 
best probe of an IGM partially ionized by X-rays prior to 
full reionization, which otherwise may not leave a unique 
imprint in the CMB. Finally, if the IGM is heated prior to 
reionization, its thermal pressure could suppress outflows 
from galaxies. A pre-heated IGM may also be detected 
by future radio telescopes in 21-cm emission against the 
CMB (Tozzi et al. 2000). 

Despite the increased electron fraction in the IGM, the 
amoimt of H^-catalyzed molecular hydrogen formed is in- 
significant, contrary to expectations. This is primarily due 
to the strong photodissociating effects of the minimum 
IR/0 and FUV (Lyman- Werner bands) backgrounds as- 
sociated with any high-z XRB on H~ and H2 respectively. 
As we have neglected radiation from stars in the treat- 
ment here, we have underestimated these photodestruction 
terms at the redshifts we consider. We therefore conclude 
that, unless the FUV and IR/0 backgrounds are strongly 
attenuated by some process which we have not included 
here (see Table 1), the positive feedback from a boosted 
electron fraction caused by a high-z XRB is not sufiicient 
to produce interesting amounts of new H2 in a uniform 
IGM. 

We mention here that there are several issues that we 
have neglected or oversimplified in this paper. First, we 
have not included the effects of any clumping in the IGM, 
which will reduce our adiabatic cooling estimates, alter the 
equilibrium calculation values of n(H2), and challenge our 
assumption that the X-rays heat the IGM homogeneously. 
Second, although we have noted the Jeans mass filtering in 
a heated IGM, we have not examined in detail the further 
cooling that must occur in virializing halos in order for 
star formation to occur. We have found that an XRB may 
not create significant amounts of new H2 in the IGM, but 
the situation may prove to be different within a halo or an 
IGM overdensity that is being irradiated by X-rays from 
a nearby source (see, e.g., Ricotti et al. (2001)). In such 
cases, the favored regimes of baryon density and tempera- 
ture for H2 production may change. It would be interesting 
to explore the feedback of a heated partially ionized IGM 
on nonequilibrium H2 chemistry in the presence of a (de- 
stroying) FUV versus (enhancing) X-ray background pho- 
ton field (Haiman et al. 2000), particularly when including 
the effects of secondary electrons in a dense, neutral halo 
at high redshift. We defer more detailed treatments of 
these problems and related consequences to future work. 
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